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I. INTRODUCTION 

Topological insulators constitute a new phase of mat- 
ter, with preserved time-reversal invariance and a bulk 
energy gap that is accompanied by topologically pro- 
tected gapless states on the boundary-^. These bound- 
ary states are helical, forming an odd number of Kramers' 
pairs where electron spin is coupled to the direction of 
propagation. On the edge of a two-dimensional (quan- 
tum spin Hall) topological insulatoi! 2 ^ 3 -!, or in a topologi- 
cal insulator nanowire^"^, such helical states will in the 
presence of electron-electron interactions give rise to the 
new class of gapless one-dimensional corre lated electron 
liquid known as the helical Luttinger liquidr^^. 

This one-dimensional electron system is distinct 
from standard spinless, spinful and chiral Luttinger 
liquids^^, and from the "spiral Luttin ger liquids" 
formed b y a sp in-orbit coupled quantum wire^ 3 -! in a mag- 
netic fiel d 14 ! 15 -! due to its time-reversal invariance with a 
single Kramers' pair. This property makes the helical 
liquid stable against weak disorder since it forbids single- 
electron backscattering from spinless impurities 1 , and al- 
though electron-electron interaction may cause impurity- 
induced correlated two-electron backscattering, the gap- 
less helical liquid is still left intact provided the screen- 
ing is not very weakPESl. However, time-reversal invari- 
ance does not protect the helical liquid against the ef- 
fects from magnetic impurities, occurring due to dopant 
ions or trapped bulk electrons caused by potential in- 
homogeneities, allowing electrons to backscatter via spin 
exchange. 

Magnetic impurities in conventional Luttinger liq- 
uids are well-studied^. Single-electron backscattering 
from non-magnetic impurities is a relevant perturbation 
causing perfect reflection for repulsive electron-electron 
interaction^. The Kondo scattering off magnetic im- 
purities leads to a low-energy behavior governed by a 
strong-coupling fixed point where the impurity is com- 
pletely screened by the electrons, approached with the 
anomalous exponents of a non-Fermi liquidMMll a nd re _ 
suiting in vanishing zero-temperature conductance. In- 
stead considering a helical Luttinger liquid, studied early 
on as a truncated model for a standard Luttinger liq- 
uid, there is only backward Kondo scattering and the 



Kondo effect was shown to map exactly onto that in a 
Fermi liquid 2 *^. This yields the typical scaling behavior 
for the Kondo effect of ordinary conduction electrons^. 
Hence, for the helical Luttinger liquid with antiferromag- 
netic Kondo coupling there is a low-temperature strong- 
coupling fixed point, where the local magnetic moment 
(stabilized against charge fluctuations by the electron 
correlations 23 ) is screened by the electr ons effectively 
leaving behind a spinless impurity^mH. The edge of 
a two-dimensional topological insulator merely bypasses 
this Kondo singlet, thus restoring unitary conductance 
G = Gq = e 2 /h at zero temperature^. Transport prop- 
erties in the vicinity of this strong-coupling fixed point 
then follow power-laws governed by the correlated two- 
electron backscattering, whereas at high temperatures 
the Kondo scattering results in a logarithmic tempera- 
ture correction in the linear conductance^, which how- 
ever vanishes in the dc limilP^. 

In addition to the atomic spin-orbit coupling that 
lies behind the very existence of the quantum spin Hall 
state^, with the idealized helical edge states with right 
(left) movers with spin up (down), there will also be ad- 
ditional spin-orbit couplings which make spin no longer 
conservecP-3 As long as the bulk gap remains, preserved 
time-rever sal invariance ensures the robustness of the he- 
lical liquid^. However, an important question is how 
the Kondo physics of magnetic impurities is affected by 
the additional spin-orbit interactions in the helical liq- 
uid. For a quantum well with structure inversion asym- 
metry (SIA), i.e. an asymmetry in the confining potential 
along the direction perpendicular to the well, there will 
be a Rashba spin-orbit interactioiPSHUl This interaction 
can be controlled by an external gate voltage, a scenario 
explored in Ref. 001 When the crystal structure lacks 
inversion symmetry, i.e. in the case of bulk inversion 
asymmetry (BI A), there will also be the Dresselhaus spin- 
orbit interactio n 129 ! 30 ! 33 -). Including this type of interaction 
hence adds a crucial piece to a fuller understanding of 
the effects of spin-orbit interactions in these systems. In 
fact, both HgTe and In As have the inversion asymmetric 
zincblende crystal structure, thereby displaying the Dres- 
selhaus interaction. HgTe/CdTe quantum wells, the first 
quantum spin Hall insulator experimentally observecP, 
is known for its large Rashba coupling^! and one here 
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expects the Dresselhaus coupling to be much smaller in 
comparison 3 -^. For InAs/GaSb quantum wells, where he- 
lical edge states also have been demo nstrate d 36 * 37 ^ , the 
Dresselhaus interaction is sizeable^ 3 ^ and realistic pre- 
dictions for spin-orbit effects must take it into account. 

It is therefore the aim of this paper to revisit the prob- 
lem of a Kondo impurity in a helical Luttinger liquid with 
spin-orbit interaction. The results announced in Ref. 0HI 
are extended to include both Rashba and Dresselhaus 
interactions, showing that the obstruction of Kondo sin- 
glet formation in certain parameter regimes persists but 
will vanish for sufficiently strong Dresselhaus interaction. 
The linear conductance far above the Kondo temperature 
is obtained, also in the low-frequency limit where a rate- 
equation approach must be employed-^. Furthermore, 
current-voltage and noise characteristics of the backseat - 
tered current are calculated, as well as the expression for 
the thermal conductance, showing how the electric-field 
adjustable Rashba interaction offers a mechanism to con- 
trol transport properties. 

II. MODEL 

The low-energy dynamics of the helical electrons, ne- 
glecting interactions, is described by the one-dimensional 
Dirac Hamiltonian 

H = v F J dx &(x) [-i<J z d x ] V(x) (1) 

where v F is the Fermi velocity, and the two-component 
spinor = (ip^,^n . Here ipf (ip±) annihilates an elec- 
tron with spin up (down) in the z direction, the growth 
direction of the quantum well. Unless stated otherwise, 
the units are such that h — ks — 1 and the lattice con- 
stant a = 1. 

Rashba and Dresselhaus spin-orbit interactions, of 
stren gth a and /3 respectively, to leading order add the 

termassnnmmi 

H R = aj dx &(x) [-ia*d x ]V(x) (2) 

H D = (3 J dx&(x)[-ia x d s ]9(x) (3) 

to the Hamiltonian. The Kondo interaction, describ- 
ing the antiferromagnetic exchange coupling between the 
electrons and a spin-1/2 magnetic impurity located at 
x = 0, is given by 

H K = Irt(O) [J x a x S x + J y a y S y + J z a z S z ] (0), (4) 

where a 1 (S l ), i = x,y, z, are the electron Pauli matrices 
(impurity spin operators). At a quantum well interface, 
spin-orbit induced magnetic anisotropy for an impurity 
would result in^21 j x = Jy ^ 

The effect of a Rashba interaction ^ on the helical 
states (fTl) was studied in Ref. H21 and that analysis can 



now be adapted to also include the Dresselhaus coupling 
Since H + Hr + Hp is given by 

J dx ¥{x) [-i {v F a z + aa v + (3a x ) d x ] *(a:), (5) 

it follows that a rotation # ->■ = e ^ v H^ e - ia " e > 2 ^ 
diagonalizes Ho + Hr + Hu — > H' , with 

H'o = v aP J dx ^(x) [-ia z d x ] Ht'{x), (6) 

where v a p = \Jvp + a 2 + /3 2 , and the rotation angles are 
determined by 

sin0 = a/^Jv F + a 2 , (7) 

sm<P = /3/^v F + a 2 + P 2 . (8) 

Thus the spectrum of Ho + Hr + Ho is linear, and the 
components and ip'_ of the spinor are left- and 
right movers. It is therefore possible to treat the effects 
of the spin-orbit interactions ([2|-([3]) exactly by working 
in the rotated basis. 

The Kondo interaction Q is now analyzed in the ro- 
tated basis. In order to have the impurity and elec- 
tron spins quantized along the same axis, rotate S — > 

S > = e <5»*/2 e -i5C0/2 Se iS-fl/2 e -iS»*/2 - Xhen the Kondo 

Hamiltonian becomes 

H' K = *' f (0) [ J x a x S' x + J' y a v S' y + J' z a z S' z 

+J xy a x S'v + J yx <J v S' x + J xz a x S' z (9) 
+J zx a z S' x + J yz o y S lz + J zy a z S' v } tf'(0), 

where the coupling constants are given by 

J' x = J x cos 2 <f> + Jy sin 2 9 sin 2 </> + J z cos 2 9 sin 2 <fi, 
J' y = J y cos 2 9 + J z sin 2 9, 

J' z = J x sin 2 4> + J y sin 2 9 cos 2 <j> + J z cos 2 9 cos 2 (j), 

Jxy — Jyx = (Jz — Jy ) cos 9 sin 9 sin <f>, 

Jx Z — J Z x — {Jx — Jy sin 2 9 — J z cos 2 9) cos <f> sin <f>, 

Jyz ~ J Z y — {Jy ~~ J z ) COS 9 Sm9 COS <f>. 

The next step is to include the electron-electron inter- 
actions allowed by time-reversal symmetry. Assuming a 
band away from half-filling, kp ^ 7r/2, Umklapp scatter- 
ing 

Km = 9um I &xe- iik ^'l{x)i,'l{x + a) 

xip'_(x)ip'_(x + a) + H.c. (10) 

can be ignorecP. Here a point splitting with the lattice 
constant has been performed. Dispersive and forward 
scattering are given by 

H' d = g d J dx ^1{xW + {x)^1{xW_{x), (11) 

(12) 
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and at the impurity site there will in a ddition be corre- 
lated two-particle backscattering^lMMI 



H 2p 



ff2p^(0)^+(a)^L(0)C(a) + H.c. (13) the Hamiltonian H 



Let us now derive the renormalization-group equations 
for the Kondo couplings in Eq. Writing the local 

partition function at the impurity site corresponding to 



HLL 



H' K as a functional integral 



and inelastic single-particle backscattering^"^ 

H'ie = <?ie(^(0)VU0)- ^(0)^.(0) 

xip1(a)ip'_(a) + H.c. 



(14) 



Employing standard bosonization method d 11 ! 12 ! the in- 
teracting Hamiltonian H' Q + H' d + H'^ forms a helical Lut- 
tinger liquid 



h 'hll 



dx [{d x ip) 2 + {d x $f 



(15) 



l/2 e -i x AF(i9± ¥ j) 



where the electron operators ip± = (2ttk 
are represented in terms of non-chiral Bose fields (p(x) 
and ~&(x), with [tp(x), d y d{y)] = iit5(x — y). Here k ss 
vp/D is the edge state penetration depth acting as short- 
distance cutoff with D the bandwidth. The velocity v = 
{(vup+gj /t) 2 — (gel/ft) 2 ] 1 ^ 2 , and the Luttinger parameter 
K = [(TTv a +g f - g d )/(irv a +g f + gd)} 1/2 - The correlated 
two-particle and inelastic single-particle backscattering 
terms in Eqs. (fT3|-(14) are now expressed as 



2(7TKj z 



H\. = 



2ir 2 VK 



: [(9^(0)] cos [V4ttK<p{0)] 



(16) 
(17) 



where : ... : denotes normal ordering. The Kondo inter- 
action ([9| is represented as 



H' 



— cos 



V[cp] exp<j- — / du \u\\<p(u)\ 2 - S K [<p(r)] 



with the imaginary time action Sk for the Kondo terms 
given by 

rP ( A 

SkMt)} = / dr \ — cos[V47rZ^(r)l 

In v 7VK 



— sin[v 4arK <p[r)\ + 



dMr)} 

(20) 



with the impurity operators A, B and C given by 
Eq. ( [19] ) . The analysis will be perturbative in the Kondo 
couplings, with So[<p] — S[tp] — Sk[<p] the unperturbed 
action. Although the J' z coupling may be treated ex- 
actly in the case of zero spin-orbit couplings^, the ad- 
ditional coupling terms introduced in the Hamiltonian 
^ by the spinor rotation precludes this line of attack. 
Since the exact analysis in J' z yields different physics 
only for J' / (ttv) > 2(K + vK), where the Kondo ef- 
fect disappears^, for reasonable strength of the electron- 
electron interaction the perturbative treatment of J' z is 
sufficient for the purposes here. 

The procedure^ is standard Wilsonian RG, where the 
field (p is divided into a slow and a fast part, ip(r) — 
ip s (r) + <pf(r), with 



=0x0(0) 



where the operators A, B and C are defined as 

A — J X S J<ntjS ^ JxzS 

B — Jy X S 



JxyS " 4 

f j' y s' v 4 

f" JzyS ' " 



u yz u 1 

^JlS' z . 



(18) 



(19) 



1 /-A/6 

Mr) = it dw (21) 



<Pf(r) 



2tt 
1 

2^ 



-A/b 



duj e-^Xw). (22) 



A/6<M<A 



An effective action if? 8 [^ s ] is then found for the slow field 
by integrating out the fast components, 



III. KONDO TEMPERATURE 



,-Ssllp.] _ 



= J T>[<Pf] 



,S[<p] 



(23) 



The low-temperature behavior of the model can now 
be determined through a perturbative renormalization- 
group analysis. First, from Eqs. (16) and (17 1 it is 



seen that the correlated two-particle and inelastic single- 
particle backscattering operators have scaling dimensions 
4K and K + 2, respectively. Hence, for Luttinger pa- 
rameter K < 2/3, the two-particle backscattering is the 
dominating of the two. For K < 1/4 it turns rele- 
vant, with a crossover from weak to strong coupling at a 



temperature^^! T 2p « Dg, 



1/(1-4K) 
2p 



Rescaling the high-energy cutoff A = v/n with the pa- 
rameter b, so that A/b is the new cutoff, gives the new 
rescaled effective action and allows extraction of the RG 
flow equations. 

A linked cluster expansion to second order in the 
Kondo couplings, 



S s [ip s ] = S [ip s ] + (S K [(p])f 

-\{{s K M 2 ) f -{s K Wf), 



(24) 
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FIG. 1: The Kondo temperature Tk as a function of the Rashba angle 6 and the ordinary Luttinger parameter Kq, for four 
different values of the Dresselhaus angle <j>. Note that the temperature scale is logarithmic, and that the vanishing Kondo 
temperature at 8, <j> = tt/2 simply reflect the diverging spin-orbit coupling strength in these limits. Left column: J x = J y > J z 
(here with J x /a — J z /a — 10 meV). Right column: J x = J y < J z (here with J x /a = 5 meV and J z /a — 50 meV). In the 
black (white) shaded area, J yz (J xz ) dominates the perturbative RG flow, hence preventing the low-temperature formation of 
a Kondo singlet. 



gives the rescaled effective actiorP^^ 
S[<p] = 5oM + 

A 



(b-b 



,2-2K 



l-K\ 



H sink/ 4itK(p(t)] 

~K 

[A,B] 



TTK 

B 

TTK 



2 V 2 TT 3 / 2 

■ \B,C] 



KVTT* 



d T ip(r) 



[VAttKip(t)] 



(25) 



+ (6-6 



1-K\ 



kvtt z ) 



where only the most relevant operators are kept. 
Among the less relevant operators that are neglected 
is the correlated two-particle backscattering operator 
cos[v / 167r Kip{r)]. It is hence generated during the RG 
flow, which is exp ected in the presence of electron- 
electron interaction d 44 l 51 [ In order to determine the 
Kondo temperature of the strong-coupling crossover, only 



the Kondo scattering which is the most relevant process 
needs to be considered. 



The commutators of the impurity operators ( 19 1 are 



[A, B] — 2% {J X yJ z ~ J xzJ zy)S X 

HJxzJzx - J' x J'z)S V (26) 
H~~ ( Jx J z y JxyJ zx)S 



[B,C] = 2i[(J y J' z -J yz J zy )S x 

~^~{JyzJzx JyxJz)^ (^'O 
~\~{^JyxJ zy ~ JyJ zx^jS 



[.A, C] — 2i (J xy J z J xz J Z y)S 

+(J XZ J ZX - J' x J' z )S y (28) 

~^~^J xJ zy J xyJ zx}$ 



Collecting the terms in Eq. (25), using Eqs. (26)- (28), 
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FIG. 2: The "phase diagram", showing which Kondo coupling 
constant that dominates the low-temperature behavior when 
Jx = Jy < Jt. [here with the same parameter values as for 
the right column in Fig. qlb] . Dark and light blue (J' x and J y , 
respectively) means the formation of a singlet between the 
electrons and impurity below the Kondo temperature. Red 
and orange (J yz and J xz , respectively) denotes the "dome" 
within which the singlet formation is obstructed since instead 
the J yz a y S' z or J xz a x S' z interactions will dominate. 



where g — gf — g<i refers to the original basis. The plots 
show a Kondo temperature that for given Kq and Dres- 
sclhaus angle <f> depends on the Rashba angle, which is 
controllable via an electric field. For "easy-plane" and 
isotropic Kondo coupling (J Xl J y > J z and J x = J y = J z , 
respectively, which lead to very similar plots) the system 
flows towards the Kondo fixed point with divergent J' x or 
J' y couplings resulting in a singlet between the impurity 
and electrons. However, for "easy-axis" Kondo coupling 
{ Jx, Jy < Jz), there are parameter regions where instead 
the couplings J yz a y S' z or J xz (r x S' z dominate the RG 
flow, preventing singlet formation. These regions are in- 
dicated in Fig. [I] but note that their precise sizes depend 
on the numerical values of J X ,J V and J 2 . In Fig. ([2| they 
are shown for all values of <fi, 6 and K , indicating that 
the "dome" is not completely symmetric under rotation 
in the rry-plane. 



IV. LINEAR CONDUCTANCE 
A. Strong-coupling limit at zero temperature 



the RG equations for the Kondo couplings become 



The low-temperature limit of the conductance is gov- 



diJ'x 


= 0--K)JL-\ 


VK ( j y J z JyzJzy) j 


Olj'y 


= {l-K)J' y -{ 


-VK{J' X J' Z - JxzJzx), 


diJ'z 


= vK{j'J y - 


JxyJyx") i 




dlJxy 


= (l-K)Jxy 


+ VK (JyzJzX ~ 


JzJyx) i 


deJyx 


= (l-K)Jyx 


+ vK(Jx Z Jzy ~ 


J z Jxy) ? 


deJxz 


= (l-K)Jxz 


-\- vJ\. i^JyxJ zy 


7' 7 ) 

<Jy<J ZX ) 


deJzx 


JjK {^JxyJyz 


- J' J ) 

• J y°XZ ) I 




diJyz 


= (±-K)Jyz 


~t~ 1/K { JxyJzx 


J xJzy) -i 


deJ Z y 


uK {JxzJyx 


— J x Jyz) j 





(29) 



in terms of the "renormalization length" £ = In 6, and 
where v = l/(nv) is the density of states. Note that 
although the RG flows are different, the bare values of 
the coupling constants, given by Eq. ^ where J xy (£ — 
0) = Jyx(£ = 0) etc., are not all independent. 



The RG equations ( 29 ) are now solved numerically. As 
the temperature T = Djb = D e~ e is reduced, the Kondo 
temperature Tk is reached when one of the Kondo cou- 
plings in Eq. |9]) becomes of the same order of magnitude 
as the energy scale of the unperturbed Hamiltonian (J6|, 
i.e. v a p. Fig. [T] shows Tr- as a function of the three in- 
dependent parameters, the Rashba angle 9, the Dressel- 
haus angle <f>, and the the ordinary Luttinger parameter 
K = K(6 = Q,<f> = 0), 



K = 



TTVp 

irvF+2g 



erned by the RG flows of the Kondo ( 18 ), correlated two- 



(30) 



particle (16) and inelastic single particle (17) backscat- 
tering terms in the Hamiltonian. The inelastic single- 
particle backscattering term H' ie has scaling dimension 
2 + K, and is hence irrelevant in the sense of RG, con- 
tributing a term SG le - T 2+2K to the conductance^ 21 . 
The correlated two-particle backscattering term H' 2p has 
scaling dimension 4K and is an irrelevant perturbation 
when K > 1/4, in which case it contributes a term 
SG bs - T SK ~ 2 to the conductance^. For K < 1/4, the 
H' 2p turns relevant, hence opening up a gap at the energy 
scale of Tb s , leading to zero conductance of the system at 
T = which is approached asP G - j"2(i/4K-i)_ 

Since the Kondo backscattering is a relevant 
(marginally relevant) perturbation for K < 1 {K = 
1), the transport properties of the system in the zero- 
temperature limit when 1/4 < K < 1 depend on what 
type of Kondo fixed point is reached. Outside the domes 
of dominating J yz a v S' z or J XZ <J X S' Z interactions, the 
impurity is locked into a Kondo singlet below the en- 
ergy scale of the Kondo temperature Tr-. For a he- 
lical edge liquid of a two-dimensional topological insu- 
lator, this means that the Kondo impurity effectively 
disappears from the problem, with no contribution to 
the zero-temperature conductance^, whereas a quantum 
wire would effectively be cut into two halves resulting 
in an insulating state. Inside the dome, where Kondo 
screening is obstructed resulting in a non-vanishing ef- 
fective spin, the diverging scattering strength suggests 
the destruction of the gapless state below Tk- 
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B. Linear response at weak coupling 

1. The ac conductance 

For temperatures T > T^, T 2p , the system is in a weak- 
coupling regime where the Kondo couplings J and the 
two-particle backscattering coupling g 2p remain small pa- 
rameters. This allows a perturbative calculation of the 
linear conductance. Since the components and 
of the spinor \J/' in the rotated basis are left- and right 
movers, the current operator I is given by 



1 = (e/2)fit (tf'Vtf') 



(31) 



Hence the part 81 of the current operator that is 
due to the Kondo scattering H' K is given by 81 = 
e/(i2K)\^'^(j z ' i S' ,H' K \. It is here possible to conveniently 
consider the effects of the J' z coupling constant in Eq. |9| 
to all orders in perturbation theory by performing the 
unitary transformation U(H' HLL +H' K )U\ withal 



(32) 



which cancels the J' z term in H' K against one that is gen- 
erated in H' HLL provided one chooses A = J^j-KvyK. 
This transforms the Kondo Hamiltonian H' K to 



H' K = 



2nn 



ijy 



ij 



2 

Jx. 



;< , . J xy 
+ I — 



'■• '- \ iV47r 



>y.v 



i^Tr{2\/~K+\)<p 



'■kK 



d^e-^^f : S' + + H.c. 



in terms of the spin raising and lowering operators S r± = 
S' x ±iS' v . The current correction operator SI then takes 
the form 



SI 



le 

2lTK 



J T + J,. 



JL - JL 



i J, 



J. 



\ i\/ Air Kip giz 



H.c. 



(34) 



Note that both H' K and 81 are local at x = 0, and for 
brevity only the time arguments of the fields are spec- 
ified. The Kubo formula for the conductance, G(u>) = 
(1/hu) J °° d* e iut ( [jt (<) , 1(0)] ) , now gives the linear con- 
ductance correction 8G = G — Go due to the Kondo scat- 
tering. To second order in the Kondo couplings (except 



for J' z which is treated exactly) , the unperturbed retarded 
Green's functions^ 



71" AC 



2 if, 



2i sin(irKj) 
( sinh |)- 2 ^ (35) 



contribute three terms in 5G 
2e 2 \A 

3 



irT 

' h 3 (2ttk) 2 [~D 



2 A', 



sin(7r Kj 



1 



dt 







(sinh^p) 2 ^ 



(36) 



where A 1 = {J' x +J' y )/2, A 2 = (J' x -J' y )/2+i{J xy +J yx )/2, 



A 3 = (iJ yz + J xz )/2, K 1>2 = K T = (VK T A/2) 2 and 
K3 = K. As in Refs. l24l and \2~E\ the standard integrals 
arc calculated by introducing a cutoff, then changing in- 
tegration variable to x — irTt/h and then to u — tanhx, 
giving in the limit u) -C T 



1 



dt 



rrTt\2Ki 



(sinh^l) 





h 1 r°° dx e-K**/-*?)* - 1 
ttT iw J ($mhx) 2K i 

2 poo 



h 
n 

ttT 



dx 



x 





2 rl 



(sinhx) 



2A, 



du 



-2K, 



(1 - M 2 ) A ' J_1 arctanhu, (37) 







which is evaluated to give 

.2 3 



(33) SG(u) = 



e 
2T 



2Kj-l 



p[ {ttFi) 2 vk T(2K j ) \ D J 



(38) 

Hence 8G(uj) is ^-independent in the limit J 2 <C w <C T, 
where J denotes the set of Kondo couplings. Rewriting 
Eq. p8| gives the final expression for the conductance 
correction from the Kondo scattering when J 2 <C w <C T 



8G = - 



1 



J x + J y\ [r(i^T-)] 2 (2nT\ 
T(2K_) \ D J 



2K- 



1/ . J xy 
+ I — 



>yx 



[T{K+)f S2kT\ 2K 
T{2K + ) { D J 

2K-2 



[T(K)} 2 ( 2-kT 
Y{2K) \~D~ 



2. The dc limit 



(39) 



It is important to keep in mind that Eq. (39) for the 



conductance correction was derived perturbatively in the 
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Kondo couplings J. Hence when lj <C J 2 the above type 
of analysis is no longer valid, as noted in Ref. [26] where 
instead a rate equation was set up for the impurity spin. 
In Ref. HOI this rate-equation approach was extended to 
incorporate the additional scattering types present when 
a Rashba interaction is added. Since a Dresselhaus inter- 
action merely changes the different rates for the scatter- 
ing events in these rate equations, the calculations in the 
supplemental material of Ref. HHI are easily modified to 
include also Dresselhaus interaction. For completeness, 
let us nevertheless carry them out explicitly also here. 

Denoting the rates for the scattering events encoded in 
the Hamiltonian Q with y± for a T S /zt , y' T for a ± S /± , 
y± for a T S' z and 75 for a z S' ± , the rate equation for the 
impurity spin becomes 

d t P+ = ( 7+ + 7- + %)P- - (7- + l'+ + %)P+, (40) 

with P± the probability of the impurity spin being in the 
+ or — state, and the probabilities obeying P + + P_ = 1. 
When adding a voltage bias V = Vbe _lw * between the 
left- and right-moving electrons, encoded by the term 



H v = 



eV 



(41) 



in the Hamiltonian, the voltage-dependent rates are given 

by 7± = 7oA±,7± = 7o A ±>7± = 7o A ±> witn A ± 
1 ± eV/2T for eV -C T. The current correction from the 
impurity scattering can now be written as 



SI = -e ( 7+ P_ - -/LP- + i+P+ - 1-P+ 
+7-/2-7^/2). 



(42) 



Substituting the solution of Eq. (40) into Eq. H2L 



one finds the linear conductance correction 5G(ufj = 
SI/(V e~ iuit ) to be 

e 2 

5G(u)) = -2y[(7o+7o + 7o)^ + *87o7o 

+i2{ lQ + 7o )( 7o z + Yo) + «27o7o] 
x[i2(7 + 7^+7o)+H _1 . (43) 

In the dc limit, i.e. when w « J 2 « T, this becomes 



SG(u -> 0) 



2 r 



e 
2T 



47o7o + (7o +7 ) (7o +7o)+7o7o 



7o+7 +7o 



(44) 



The remaining task is to determine the rates 70, 7 , 7q 
and 7n ■ In the high-frequency regime one can compare 
Eq. {43} , which gives 8G{lo > 7) = -{e 2 /2T)(jo + 7o + 



7q), with the linear-response result in Eq. (39 1. This 
gives 



To 

7o 
7o Z 



(4 



^{r/DfWK-xnr-i 

\j' x -j' v + i{J x y + J VX )\\T/Df^+W 2 
\lJy Z + J XZ \ 2 (T/D) 2K - 1 



(45) 



In order to obtain the rate 75* of impurity spin-flips not 
accompanied by electron backscattering, one needs to do 
an analogous linear-response analysis for a field that in- 
stead couples to the impurit}^]. Such a calculation yields 



\tJ zy + J zx \ 2 T/D 



coming from the last term in Eq. (133 



(46) 



V. BACKSCATTERING CURRENT 
A. I-V characteristics 

A finite time-independent bias voltage V = Vq is con- 
veniently treated in a picture where the voltage term (41 ) 



in the Hamiltonian is replaced by a time dependence in 
the electron fields, tp± — > e~ iHvt ^ ti il)L. In this picture 



H' K (t) 



1 



2lTK 

+A 2 e ieVt / h e i Vn(2*/K+\)<p gf- 

+A 4 -L= : d x de' l ^ Xip : S' + 
VttK 



(47) 



+ H.c. 



and 



SI(t) = 



2nn 



JeVt/H i^(2y/K+\)ip at- 



(48) 



+A3 e ieVt / h e * -S^wKip g/z 



H.c. 



where Ax = (J' x +J y )/2, A 2 = (J' x -J> y )/2+i{J xv +J yx )/2 
and A3 = (iJ yz + J xz )/2 as before, and A4 = (iJ zy + 
J zx )/2. 

The non-equilibrium expectation value of the current 
operator is given by 

(^W) = \ J2( f K{si(^)e- l ^ dt ' H '-^}), (49) 

r)=± 

where Tk is the time ordering operator on the Keldysh 
contour K, and t is time t either on the upper (+) or 
lower (— ) branch of K, see e.g. Ref. 1551 To lowest order 
in the Kondo couplings, the most relevant terms are 



(6I(t)) = 



(S'+S'-) 



x ( ^(fK^V^^^^e^v^ 1 ^^ + H.c.}) 

+ |A2| 2 (Tx{e l \ /3 ^+ v(t ' ,) e-V 4 ^+^(*') + H.c.}) ) (50) 

+ \A 3 \ 2 (f K {e lV ^ v( - t ^e- lV ^ v(t '^ + H.c.})(S' z S' z ) 




FIG. 3: The backscattered current SI due to the Kondo impurity, for different values of the Rashba and Dresselhaus angles 
9 and <j> (parameterizing the strengths of the spin-orbit interactions), and the Luttinger parameter Ko- Left: Current- voltage 
characteristics for 81 at fixed temperature T. Right: 81 at fixed voltage and temperature, plotted as a function of the Rashba 
angle 6. The coupling constants have been RG improved at the fixed temperature. With typical parameter values as discussed 
in the text, voltages on the order of 1 mV at temperatures in the mK regime will give current corrections 81 on the order of a 
nA. 



The integrals on the Keldysh contour are founcP^ to re- 
sult in 



(81®) 



2tt 2 vk 



eV 
2T 



sinh i ]^ 



\T{K 3 +ieV/2T:T) 
T(2KA 



\A,\ 2 

2 ' 



2ttT\ 



2if,--l 



(51) 



where T is the gamma function. This could equivalently 
be written in terms of the Euler beta function B as 



51 



2ttT 



2Kj-l 



2ir 2 vK^~"l V D 

xB(Kj + ieV/2irT, Kj - ieV/2irT) 
smfrjKj - ieV/2nT)] X 



cos(n Kj) 



which is on the form of the result p resent ed in Ref. 1401 
derived using the tunneling formalis m 1 54 ! 55 !. Note the dif- 
ference in sign of the prefactors. Since the backscattering 
current SI = I — Iq, where Iq is the current in the ab- 
sence of the impurity, it is negative. The current correc- 
tions due to correlated two-particle, and inelastic single- 
particle, backscattering can be found in Ref. |46] 

In Fig. [3] the backscattered current SI, given by 
Eq. ( |5"T| ) , is plotted for different strengths of the Rashba 
and Dresselhaus interactions and Luttinger parameters 
Kq. As an illustration, possible parameter values can 
be collected from Ref. HO] for the case of an Mn 2+ ion 
implanted in the helical edge liquid of a HgTe quan- 
tum welPR For this material,^ vf ~ 5.0 x 10 5 m/s, 
the lattice constant 57 a ps 0.5 nm, and the bandwidth 3 
D m 300 meV. Estimates for the value of the Lut- 
tinger parameter Kq for the electron-electron interaction 



strength vary betweeiP^ESHU] q.5 and 1, depending on 
the geometry and composition of the specific semicon- 
ductor heterostructure. The Kondo coupling constants 
are taken to be J x = J y = 3J Z — 3 J/, where a rough 
estimate yields^HUMl j I j a on the order of 10 meV. For 
HgTe quantum wells, values of the Rashba coupling ha 
are typically in the range betweerPS 5 x 10 -11 eVm to 
1 x 10 -10 eVm, i.e. the value of the Rashba rotation angle 
is roughly in the interval 0.14 <#< 0.27 for typical sam- 
ples and voltages. For this particular heterostructure, the 
Dresselhaus interaction strength is typically negligible 
compared to the Rashba 3 - 5 -!, but note there are other ma- 
terials with helical edge states where Dresselhaus is com- 
parable in mag nitude or even stronger than the Rashba 
interactiorPSES. With these parameter values, a voltage 
bias in the mV regime at a temperature around 1 mK 
typically yields a current SI on the order of 1 nA, i.e. a 
correction of roughly one percent to Iq = GqV . 



B. Current fluctuations 

From the above results it is an easy task to extract 
the noise in the current correction SI. In terms of the 
symmetric combination of the noise correlators 

S(t-t') = (si{t)si{t')) + (si(t')si{t)) 

-2(Si(t'))(SI(t)), (53) 
which in the Keldysh formalism takes the form 



S(t-t') = ^2{f K {<5/(t")57(t'-")e- i ^ dt ' K « (t ')}) 

7J=± 



-2(SI(t)), 



(54) 



the noise spectrum S(u) follows as the Fourier trans- 
form S(oj) = dse MS S(s); for a review see Ref. [S3] 
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To lowest order in the Kondo couplings, Eq. (54) 
reduces to a current-current correlator, S(t — t') — 
E n =±<2V{W(i 7? )^(i'"' ? )})-2(W(i)), which is evaluated 



precisely as Eq. ( 50 ). This yields the zero-frequency noise 



S(lu = 0) = 



cosh 



eV 
2T 



3=1 L 



2Kj 



|r (Kj + ieV/2nT) 
WW) 



(55) 



i.e. S(u> = 0) = 2e\(5I(t))\coth(eV/2T). Since 
coth(el^/2T) sa 1 for the temperature and voltages used 
in Fig.[3j this gives a simple relation between the current 
and its zero-frequency noise. 



VI. THERMAL TRANSPORT 

Let us also consider the thermal conductance of the 
system. In terms of left-and right-moving fields 4>l,r = 
(i? ± ip)/2, such that H' HLL = (hv/2) J dx[(d x cj) R ) 2 + 
{dx^h) 2 } = H'hll + H'hll> tne backscattered thermal 
current SIq is given by the operato:^ 



IV 

~2 



which gives 



SIq = x Lp8l 



&xd t {H'g LL -H'£ LL )/2 (56) 
dx[:d x (pd x #:(x), H' K (0)}, 



(57) 



+A 4 ^1 ( : flgtfe-V^ : S'+ + H.c.) 
2^K V * J 



where H' K and SI are given by Eqs. (33) and (34), and 



Aj = (iJ zy + J zx )/2. The first term in Eq. (57) has 



its analog in Refs. [M] and [251 but the second term is 
new and is a result of the spin-orbit interaction. The 
correction SIC to the linear thermal conductance due to 
the Kondo impurity is obtained from the Kubo formula, 
SIC = (1/hw) J °° dte iut {[Sp Q (t),8i Q (0)]), where the two 
terms in Sin add separately to 5)C. Now, the contribu- 



tion from the first term in Eq. (57) follows from Ref 



The contribution from the second term follows by noting 
the close resemblance between that operator, with scal- 
ing dimension 2 + A 2 /4, and the inelastic single-particle 
backscattering term in Eq. (17), with scaling dimension 
2 + K, which makes it possible to use known resultiP^. 
The two contributions sum to give the linear thermal 



conductance correction 



8e 2 k B T 2 



did 



uj 2 ReSG{uj) 
smh 2 {hu: / 2k bT) 

2ttT\ x2/2+3 



■ 



AK K 3 k B T 2 V D 
'(l + 22A 2 /4)(A 2 /4 
(A 2 /2 + 3)(A 2 /2 + 
[r(A 2 /4 + l)] 2 
T(A 2 /2) ' 



1) 



9A^ 
4 



(58) 



Comparing the temperature dependence of 5G(w) in 
Eq. (43} with the factor - (T/D) x2 / 2+3 in the last term 
in Eq. ([58]), it follows that to leading order the usual 
formula 



SIC 



h 3 K 
8e 2 k B T 2 



dio 



lu 2 Rc SG{u) 
sinh 2 (huj/2k B T) 



(59) 



is valid also here. The same observation as in Ref. [5S] 
therefore applies, namely that the thermal conductance 
SK, is determined by the electrical conductance SG(uj) 
within a frequency range < lu < 2k B T/h, and that even 
when SG vanishes in the dc limit this is not necessarily 
true for SIC. 



VII. DISCUSSION 

The interplay between a magnetic impurity and spin- 
orbit interaction in a helical Luttinger liquid has been 
studied. Using a transformatiorPS that effectively re- 
moves the spin-orbit interactions among the electrons, at 
the price of introducing additional coupling terms to the 
impurity, the results are non-perturbative in the spin- 
orbit couplings. The effect described in Ref. [301 that 
a Rashba spin-orbit interaction appear to prevent the 
Kondo singlet formation in certain parameter regimes, is 
found to persist also with a Dresselhaus spin-orbit inter- 
action. The perturbative RG analysis therefore suggests 
that a time-reversal invariant perturbation can destroy 
the Kondo effect in one dimension, in stark contrast to 
the c ase fo r ordinary metals where the Kondo effect is 
robust 65 and such perturbations only shift the Kondo 
temperature^^. 

In the other parameter regimes, the Kondo effect will 
set in. Since the Rashba interaction can be controlled 
by an external gate voltage, this offers a mechanism to 
electrically control the Kondo effect and therefore the 
transport properties of the edge at low temperatures^. 
At higher temperatures the Rashba dependence of the 
Kondo scattering permits a way to electrically control the 
small current corrections. For in-plane symmetric Kondo 
couplings J x = Jy 7^ J z , and negligible Dresselhaus in- 
teraction, Eq. (44) reveals a vanishing dc linear conduc- 



tance correction which can be made non-zero by tuning 
the gate voltage. With non-zero Dresselhaus interaction 
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strength, the linear conductance correction however re- 
mains non-zero also in the dc limit even without Rashba 
interaction. 

Unfortunately, there is no exactly solvable Toulouse 
limit as in the absence of spin-orbit interactions^. This 
rules out an exact analysis of the low-frequency behavior 
at some special value for the interaction strengths, and 
has its root in the new operators generated by the spin- 
orbit interaction bringing several different scaling dimen- 
sions into the problem. The different exponents show up 
in the temperature dependence of conductances, and in 
the temperature and voltage dependence of the current 
correction due to the Kondo impurity. 

The spin-orbit interactions hence introduce new types 
of scattering mechanisms compared to those for the orig- 
inal helical liquid with a Kondo impurity, adding yet an- 
other piece to the different types of correlation effects 
appearing in topological insulator materials^, including 



Kondo screening j 9 l 24 H 2£l, trans port propertie d 24 * 26 * 64 ! and 
current noise correlation a 73 * 74 ! There is also the possibil- 
ity of exotic multi-channel Kondo physics, like the two- 
channel effect appearing with a quantum dot between 
two helical edges-^ where the effects of a Rashba cou- 
pling were studied very recently^. Interesting extensions 
include the case of many impurities^ forming a Kondo 
latticed on a helical edge, where the effects of spin-orbit 
interactions remain unexplored. 
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